Gap and subgap tunnelling in cuprates 
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We describe strongly attractive carriers in cuprates in 
the framework of a simple quasi-one dimensional Hamiltonian 
with a local attraction. In contrast with the conventional BCS 
theory there are two energy scales, a temperature independent 
incoherent gap Ap and a temperature dependent coherent gap 
Ac(r) combining into one temperature dependent global gap 
A = (Ap + Ac)^^^. The temperature dependence of the gap 
and single particle (Giaver) tunnelling spectra in cuprates are 
quantitatively described. A framework for understanding of 
two distinct energy scales observed in Giaver tunnelling and 
electron-hole reflection experiments is provided. 

PACS numbers:74.20.-z, 74.65. +n,74.60.Mj 



There is convincing experimental evidence that the 
pairing of carriers takes place well above Tc in under- 
doped cuprates (for a review see Ref. [0). If carriers are 
paired their magnetic moments compensate each other so 
one could expect that the normal state uniform magne- 
tization should fall with decreasing temperature because 
more and more holes are bound into singlet pairs. This 
unexpected drop of the normal state magnetic suscepti- 
bility was experimentally observed and explained in 
the framework of the bipolaron theory of cuprates . 
There is also a gap in the tunnelling and photoemission, 
which is almost temperature independent below Tc |^] 
and exists well above Tc , so that some segments of 
a 'large Fermi surface' are actually missing l9[[l0[|. Ki- 
netic [|ll| and thermodynamic |l^] data suggest that the 
gap opens in both charge and spin channels and exists 
at any relevant temperature in a wide range of doping. 
A plausible explanation is that the normal gap is half 
of the bipolaron binding energy although alternative 
explanations have also been proposed. The temperature 
and doping dependence of the gap still remains a sub- 
ject of controversy. Moreover, reflection experiments, 
in which an incoming electron from the normal side of 
a normal/superconducting contact is reflected as a hole 
along the same trajectory [|l3|, revealed a much smaller 
gap edge than the bias at the tunnelling conductance 
maxima in a few underdoped cuprates |l4| . Recent in- 
trinsic tunnelling measurements on a series of Bi '2212' 
single crystals |15|| showed distinctly different behaviour 
of the superconducting and normal state gaps with the 
magnetic field. Such coexistance of two distinct gaps in 
cuprates is not well understood [ p^p^ . 

In this letter we propose a model, which describes the 
temperature dependence of the gap, tunnelling spectra 
and electron-hole reflection in cuprates. The assump- 
tion is that the attraction potential in cuprates is large 



compared with the Fermi energy. The main point of our 
letter is independent of the microscopic nature of the at- 
traction. Real-space pairs might be lattice and/or spin 
bipolarons 0|, or any other preformed pairs. 

We start with a generic one-dimensional Hamiltonian 
including the kinetic energy of carriers in the effective 
mass (m) approximation and a local attraction potential, 
V{x - x') = -US{x - x') as 
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where s =t, i is the spin [fi = fc^ = 1). The flrst band 
to be doped in cuprates is the oxygen band inside the 
Hubbard gap as established in polarised photoemission 
|l7| , p^ . This band is almost one dimensional as discussed 
in Ref. , so that our (quasi) one-dimensional approx- 
imation is a realistic starting point. 

Solving a two-particle problem with the (5-function po- 
tential one obtains a bound state with the binding energy 



2A. 



1 
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and with the radius of the bound state r = 2/mU . We 
assume that this radius is less than the inter-carrier dis- 
tance in cuprates, which puts a constraint on the dop- 
ing level, Ep < 2Ap, where Ep is the free-carrier Fermi 
energy. Then real-space pairs are formed. If three- 
dimensional corrections to the energy spectrum of pairs 
are taken into account (see, for example, Ref. [^) the 
ground state of the system is the Bose-Einstein conden- 
sate. The chemical potential is pinned below the band 
edge by about Ap both in the superconducting and nor- 
mal state [Ql , so that the normal state single-particle gap 
is Ap. The binding energy 2Ap might change due to the 
same corrections. However, this change does not affect 
our further results as soon as they are expressed in terms 
of Ap rather than U. 

Now we take into account that in the superconduct- 
ing state (T < Tc) the single-particle excitations interact 
with the condensate via the same potential U. Applying 
the Bogoliubov approximation |2l] we reduce the Hamil- 
tonian, Eq.(l) to a quadratic form as 
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dx[Aci^\{x)i)l{x) + H.c], 
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where the coherent pairing potential 

^c = ~U{^l{x)i;^{x)) 



(4) 



is proportional to the square root of the condensate den- 
sity, Ac = constant x rio(T)^/^. The single-particle exci- 
tation energy spectrum E{k) is found using the Bogoli- 
ubov transformation as 



E{k) = [(fcV2TO+Ap)2-f A2] 
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if one assumes that the condensate density does not de- 
pend on position. This spectrum is quite different from 
the BCS quasiparticles because the chemical potential is 
negative, /i = — Ap. The single particle gap. A, defined 
as the minimum of E{k), is given by 



A = [A^ + A^] 
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It varies with temperature from 

1 /2 

A(0) — [Ap + Ac(0)^] at zero temperature down to 
the temperature independent Ap above Tc. The conden- 
sate density depends on temperature as 1 — (T/TcY^^ 
in the ideal three (d = 3) and (quasi) two-dimensional 
{d = 2) Bose-gas. In the three-dimensional charged 
Bose-gas it has an exponential temperature dependence 
at low temperatures due to a plasma gap in the Bogoli- 
ubov collective excitation spectrum [p2[, which might be 
highly anisotropic in cuprates Near Tc one can ex- 
pect a power law dependence, nQ{T) oc 1 — (T/Tc)" with 
n > d/2 because the condensate plasmon Q depends on 
temperature. The theoretical temperature dependence, 
Eq.(6) describes well the pioneering experimental obser- 
vation of the anomalous gap in Y Ba2Cu307^s in the 
electron-energy-loss spectra by Demuth et al 
with Ac(r)2 = Ac(0)2 X [1 - (T/Tc)"] below 
above Tc, and n — A. 

The normal metal-superconductor (SIN) tunnelling 
conductance via a dielectric contact, dl/dV is propor- 
tional to the density of states, p{E) of the spectrum 
Eq. (5) . Taking also into account the scattering of single- 
particle excitations by a random potential, thermal lat- 
tice and spin fluctuations one finds at T = |19 



2|, Fig.l, 
and zero 



di/dV = constant x [p 
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A is the asymmetry coefficient |T^, Ai{x),Bi{x) the 
Airy functions, and eo is the scattering rate. We com- 
pare the conductance, Eq.(7) with one of the best STM 
spectra measured in A^i-substituted Bi2Sr2CaCu20s+x 
single crystals by Hancottee et al ||], Fig. 2. This ex- 
periment showed anomalously large 2A/Tc > 12 with 



the temperature dependence of the gap similar to that 
in Fig.l. The theoretical conductance, Eq.(7) describes 
well the anomalous gap/Tc ratio, injection/emission as- 
symmetry, zero-bias conductance at zero temperature, 
and the spectral shape inside and outside the gap region. 
There is no doubt that the gap. Fig. 2 is s-like, which is 
compatible with the phase-sensitive experiments in 
the framework of the bipolaron theory [|l9|. Within the 
theory the single-particle gap might be almost k indepen- 
dent while the symmetry of the Bose-Einstein condensate 
wave-function (i.e. of the order parameter) is d — wave. 

Finally, we propose a simple theory of the tunnelling 
into bosonic (bipolaronic) superconductor in the metal- 
lic (no-barrier) regime. As in the canonical BCS ap- 
proach applied to the normal metal-superconductor tun- 
nelling by Blonder, Tinkham and Klapwijk p5| and to 
the normal-superconductor boundary in the intermedi- 
ate type I state by one of us ||l^, the incoming elec- 
tron produces only outgoing particles in the supercon- 
ductor {x > I), allowing for a reflected electron and (An- 
dreev) hole in the normal metal (x < 0). There is also 
a buffer layer of the thickness I at the normal metal- 
superconductor boundary ( a; = 0), where the chemical 
potential with respect to the bottom of the conduction 
band changes gradually from a positive large value fi in 
the metal to a negative value — Ap in the bosonic super- 
conductor. We approximate this buffer layer by a layer 
with a constant chemical potential fib (— Ap < ph < p) 
and with the same strength of the pairing potential Ac as 
in the bulk superconductor. The Bogoliubov-de Gennes 
equations may be written as usual [ p5[ , with the only 
difference that the chemical potential with respect to the 
bottom of the band is a function of the coordinate x, 
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Thus the two-componet wave function in the normal 
metal is given by 

< 0) - J) e*«^^ ^ + " (l) 

(10) 

while in the buffer layer it has the form 

^-6(0 <x <l)^a( e'P^'^ aI 

+ 7 ( i ) e-*f^" + 5(^1^^ c'P'\ (11) 

where the momenta associated with the energy E are 

Q± = [2m(/i±S)]i/2 (12) 

and 

p^ = [2m{p,±0Y'^ (13) 
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with ^ = (£'2 - A2)1/2_ The well-behaved solution in the 
superconductor with negative chemical potential is given 
by 

>l) = c(^^^ e**-'^" + d(^_A^^ e*'="^ (14) 

where the momenta associated with the energy E are 

fc± = [2m(-Ap±e)]'/'. (15) 

The coefficients a,b, c, d, a, P,"f, 6 are determined from 
the boundary conditions, which are continuity of ■ip{x) 
and its derivatives at a; = and x = I. Applying the 
boundary conditions, and carrying out an algebraic re- 
duction, we find 

a-2A,g+(p+r.9+-p-/+.g-)/^, (16) 



localised single-particle impurity states and into incoher- 
ent ('supracondensate') bound pair states, which might 
explain a significant featureless background in the subgap 
region [Q. 

In conclusion, we have proposed a simple general 
model, which provides an explanation of the tempera- 
ture dependence of the gap and of the single-particle tun- 
nelling spectra in cuprates. The main assumption is that 
the attractive potential is large compared with the Fermi 
energy, so that the ground state is the Bose-Einstein con- 
densate of tightly bound pairs. We have developed a the- 
ory of tunnelling in the metallic regime with no barrier 
and found two different energy scales in the transmission 
as observed experimentally. 

We acknowledge support of this work by the Lever- 
hulme Trust (London), grant VP/261. 



6 = -1 + 2q+[[E + Of+{q-r -p-g-) 
-{E~0r{Q-r~P+9+)]/D, (17) 

with 

D = {E + Oiq^.f++p+g+){q-f- -p-g-) 

- {E-^){q+f- +p-g-){q-f+ -p+g+), (18) 

and/* = p^ cos{p^l) — ik^ sm{p^l), —k^cos{p^l) — 
ip^ sin(p*Z). 

The transmisson coefficient for electrical current, 1 -|- 
|ap — |5p is shown in Fig. 3 for different values of I when 
the coherent gap Ac is smaller than the pair-breaking 
gap Ap, and in Fig. 4 for the opposite case, Ap < Aj,. In 
the first case, Fig. 3, we find two distinct energy scales, 
one is Ac in the subgap region due to electron-hole reflec- 
tion and the other one is A, which is the single-particle 
band edge. On the other hand, there is only one gap 
Ac, which can be seen in the second case. Fig. 4. We 
notice that the transmission has no subgap structure if 
the buffer layer is absent {I = 0) in both cases. In the 
extreme case of a wide buffer layer, I >> {2mAp)~^/^ , 
Fig. 3, or I >> (2r7iAc)~^/2, Fig. 4, there are some oscilla- 
tions of the transmission due to the bound states inside 
the buffer layer. It was shown in Ref. Q| that the pair- 
breaking gap Ap is inverse proportional to the doping 
level. On the other hand, the coherent gap Ac scales 
with the condensate density, and therefore with the crit- 
ical temperature, determined as the Bose-Einstein con- 
densation temperature of strongly anisotropic 3D bosons 
. Therefore we expect that Ap >> Ac in the under- 
doped cuprates. Fig. 3, while Ap < Ac in the optimally 
doped cuprates. Fig. 4. Thus the model accounts for the 
two different gaps experimentally observed in Giaver tun- 
nelling and electron-hole reflection in the underdoped 
cuprates and for a single gap in the optimally doped sam- 
ples ||l^. The transmission. Fig. 3 and Fig. 4, is entirely 
due to the coherent tunnelling into the condensate and 
(or) into the single-particle band of the bosonic super- 
conductor. There is also an incoherent transmission into 
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Figure Captures 
Fig.l. Temperature dependence of the gap, Eq.(6) 
(line) compared with the experiment |^^(dots) for Ap = 
0.7A(0) . 

Fig. 2. Theoretical tunnelling conductance, Eq.(7) 
(line) compared with STM conductance in Ni-doped 
Bi2Sr2CaCu20s+x §) (dots) for 2A = 90 meV, A = 
1.05, eo = 40 meV. 

Fig. 3. Transmission versus voltage (measured in units 
of Ap/e) for A^ = 0.2Ap, n = lOAp, fib = 2Ap and Z = 
(thick line), I = 1 (thick dashed line), I = 4 (thin line), 
and I = 8 (thin dashed line) (in units of l/(2mAp)^''^). 

Fig. 4. Transmission versus voltage (measured in units 
of Ac/e) for Ap = 0.2Ac, fi = lOA^, fib = 2Ac and I = 
(thick line), I = 1 (thick dashed line), I — 4 (thin line), 
and I = 8 (thin dashed line) (in units of l/(2mAc)"^/^). 
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